Dynamical singlets and correlation-assisted Peierls transition in VO2 
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A theory of the metal-insulator transition in vanadium dioxide from the high- temperature rutile 
to the low- temperature monoclinic phase is proposed on the basis of cluster dynamical mean 
field theory, in conjunction with the density functional scheme. The interplay of strong electronic 
Coulomb interactions and structural distortions, in particular the dimerization of vanadium atoms 
in the low temperature phase, plays a crucial role. We find that VO2 is not a conventional Mott 
insulator, but that the formation of dynamical V-V singlet pairs due to strong Coulomb correlations 
is necessary to trigger the opening of a Peierls gap. 

PACS numbers: 71.27.+a, 71.30.+h, 71.15.Ap 



Vanadium dioxide (VO2) undergoes a first-order tran- 
sition from a high-temperature metallic phase to a 
low-temperature insulating phase at almost room- 
temperature (T — 340 K). The resistivity jumps by 
several orders of magnitude through this transition, and 
the crystal structure changes from rutile (R-phase) at 
high-temperature to monoclinic (so-called Mi-phase) at 
low-tempcrature. The latter is characterized by a dimer- 
ization of the vanadium atoms into pairs, as well as a 
tilting of these pairs with respect to the c-axis. 

Whether these structural changes are solely respon- 
sible for the insulating nature of the low-T phase, or 
whether correlation effects also play a role, has been a 
subject of much debate. The strong dimerization, as well 
as the fact that the insulating phase is non-magnetic 
suggests that V0 2 might be a typical case of a Peierls 
insulator. However, pioneering experimental work by 
Pouget et al. demonstrated that minute amounts of Cr- 
substitutions 0], as well as, remarkably, uniaxial stress 
applied to pure VO2 leads to a new phase (M2) in 
which only half of the V-atoms dimerize, while the other 
half forms chains of equally-spaced atoms behaving as 
spin- 1/2 Heisenberg chains. That this phase is also insu- 
lating strongly suggests that the physics of VO2 is very 
close to that of a Mott-Hubbard insulator. Zylbersztejn 
and Mott 0, Sommers and Doniach [j| and Rice et al. Q 
suggested that Coulomb repulsion indeed plays a major 
role in opening the insulating gap. 

The main qualitative aspects of the electronic structure 
of VO2 have been explained long ago by Goodenough • 
In the rutile structure (space group P^/rnnm), the V 
atoms are surrounded by O octahedra forming an edge- 
sharing chain along the c-axis. The d-levels of the V-ions 
are split into lower lying ti g states and e g states. The 
latter lie higher in energy and are therefore empty. The 
tetragonal crystal field further splits the t2 g multiple! into 
an dig state and an e g doublet (dji and ir* states, respec- 
tively in the terminology of Ref. [2|). The a\ g orbitals are 
directed along the c-axis, with good cx-bonding of the V-V 



pair along this direction. In the monoclinic phase (space 
group P2i/c), the dimerization and tilting of the V-V 
pairs results in two important effects. First, the a\ g (d\\ ) 
band is split into a lower-energy bonding combination 
and a higher-energy antibonding one. Second, the Vd-O p 
antibonding e g (it*) states are pushed higher in energy, 
due to the tilting of the pairs which increases the overlap 
of these states with O states. In Goodenough's picture, 
the single d-electron occupies the dy-bonding combina- 
tion, resulting in a (Peierls-like) band gap. 

Electronic structure calculations based on density 
functional theory within the local density approximation 
(DFT-LDA) have since provided support for this quali- 
tative description in terms of molecular orbitals (see the 
recent work of Eyert || for an extensive discussion and 
references). Molecular dynamics calculations by Wentz- 
covich et al. |IJ with variable cell shape successfully found 
the Mi structure to have the lowest total energy, with 
structural parameters in reasonable agreement with ex- 
periment. DFT-LDA calculations fail however to yield 
the opening of the band gap: the top of the bonding 
d||-band is found to overlap slightly with the bottom 
of the 7r*-band (only for a hypothetical structure with 
larger dimerizations would the band gap fully open). 
Not surprisingly, recent DFT-LDA calculations of the M2 
phase |a| also fail in producing an insulator. 

This discussion makes clear that only a theoretical 
treatment in which structural aspects as well as corre- 
lations within V-V pairs are taken into account on equal 
footing, is able to successfully decide on the underlying 
mechanism for the metal- insulator transition in VO2. In 
this Letter, we fulfill this goal by using a cluster exten- 
sion of dynamical mean-field theory (C-DMFT) ^J] i n 
combination with state-of-the-art DFT-LDA calculations 
within the rece ntly developed Nth-order muffin-tin or- 
bital (NMTO) 11] implementation. This allows for a 
consistent description of both the metallic rutile phase 
and the insulating monoclinic phase. We find that the in- 
sulating state can be viewed as a molecular solid of singlet 
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dimers in the Heitler-London (correlated) limit. While a 
description in terms of renormalised Peierls bands is pos- 
sible at low-energy, broad Hubbard bands are present at 
higher energy. Recent photocmission experiments can be 
successfully interpreted on the basis of our results, which 
also yield specific predictions for inverse-photoemission 
spectra. 

Previous theoretical work on VO2 based on DMFT has 
recently appeared 0, 0. These works, however, are 
based on a single-site DMFT approach. This is appro- 
priate in the metallic phase, but not in the insulating 
phase. While we do find that by increasing U to un- 
physically high values, a Mott insulator can be induced 
in a single-site DMFT approach, the formation of singlet 
pairs resulting from the strong dimerization can only be 
captured in a cluster extension of DMFT in which the 
dimers are taken as the key unit. Only then can a non- 
magnetic insulator with a spin-gap be obtained, in agree- 
ment with experiments. C-DMFT allows for a consistent 
extension to the solid-state of the simple Heitler-London 
picture of an isolated molecule. Electrons can be shared 
between all singlets through the self-consistent electronic 
bath, resulting in a "dynamical singlets" description. 

Given the filling of one d-electron per vanadium, and 
the crystal- field splitting separating the ti g and e" g states 
in both phases, it is appropriate to work within a set 
of localized V-centered t^ g Wannier-orbitals, and to ne- 
glect the degrees of freedom from all other bands. As in 
Ref. > our Wannier orbitals are orthonormalized NM- 
TOs, which have all partial waves other than V- d xy , d yz 
and d zx downfolded. These notations refer to the local co- 
ordinate axes (z || [110]) attached to a given V atom sur- 
rounding the oxygen octahedron. DFT-LDA calculations 
followed by this downfolding procedure yield a Hamilto- 
nian matrix H^^(k) (of size 6 x 6 in the rutile phase, 
which has 2 V-atoms per unit cell, and 12 x 12 in the mon- 
oclinic phase, with 4 V atoms per unit cell). Our results 
for the DFT-LDA electronic structure in both phases are 
in agreement with previously published studies, and with 
the qualitative picture described above. The partial den- 
sities of states for the a\ g = d xy and e g = {d yz , d xz } are 
presented in Fig. and Fig. [3 for the rutile and mon- 
oclinic phases, respectively. The total ti g bandwidth is 
almost the same for both phases (~ 2.59eV in rutile and 
~ 2.56eV in monoclinic). In the rutile phase, the sin- 
gle d-electron is almost equally distributed between all 
three orbital components within LDA (0.36 in a\ g and 
0.32 in each of the e g 's). From the Hamiltonian matrix 
in real space, we can extract the hopping integrals be- 
tween V centers. The largest one, t xy , xy = — 0.31eV is 
between a\ g orbitals forming chains along the c-axis, but 
the hoppings t xz>yz within the chains and t' yzyz ,t' xzxz 
are only twice smaller (0.17 — 0.19eV). Given that there 
are 8 next-nearest neighbours of this type, and only 2 
along the chains, the properties of the rutile metal are 
therefore fairly isotropic. The electronic structure dras- 



tically changes in the monoclinic phase (Fig. 0). The 
a\ g state is split into a bonding combination below the 
Fermi level, and an antibonding combination higher in 
energy. The splitting between the bonding and antibond- 
ing states is set by the intra-dimer hopping, which we 
find to be t^xy = — 0.68 eV, hence a splitting of order 
2txy,xy ~ 1-4 eV. The two a% g peaks in the density of 
states (DOS) are very narrow, corresponding to the very 
small inter-dimer hopping f X y° xy ~ —0.03 eV. The intra- 
dimer hopping is by far the dominant one in this phase, 
with the second largest being t' xz yz in the [111] direc- 
tion (~ 0.22 eV) and all others much smaller. Also, the 
e g states are pushed higher in energy in the monoclinic 
phase, so that the LDA occupancies are now 0.74 for the 
ai g (bonding) band and only 0.12 and 0.14 for the d yz 
and d xz , respectively. Still, these bands overlap weakly 
and the LDA fails to open the insulating gap (~ 0.6 eV 
experimentally) . 

We use the LDA-NMTO Hamiltonian ff^(k) as a 
starting point for the construction of a multi-band Hub- 
bard Hamiltonian of the form of Eq. (1) in Ref. 0] 
involving direct and exchange terms of the screened 
on-site Coulomb interaction U mm > and J mm ', with the 
parametrization U mm — U, U mm ' = U — 2 J and J m m' — 
J for m^m! . We assume double counting corrections 
to be orbital-independent within the t-z g manifold, thus 
resulting in a simple shift of the chemical potential. Re- 
cently, it has become feasible to solve this many-body 
Hamiltonian using the dynamical mean-field approxima- 
tion (DMFT) and to obtain realistic physical properties, 
even when all off-diagonal terms in orbital space in the 
local self-energy £ mm ' are retained [l4|. Here, we go 
one step further by including also non-local terms in the 
self-energy. The latter are constructed from a cluster 
LDA+DMFT treatment More precisely, instead of 
calculating the self-energy from a local impurity model 
embedding one single atom in a self-consistent bath, a 
pair of V atoms in the complex bath is explicitly con- 
sidered. The self-energy £^j; and (Weiss) dynamical 
mean-field G^m' become matrices in both the orbital in- 
dices m, ml — (xy, yz, xz) and the intra-dimer site indices 
icjc — 1, 2. The inter-dimer components of the self- 
energy as well as long-range correlations are neglected in 
this C-DMFT scheme. Using the crystal symmetries in 
the monoclinic phase, a 12 x 12 block-diagonal self-energy 
matrix is constructed from the 6x6 matrix £ so that 
the C-DMFT approximation to the self-energy takes in 
this case the form : 
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where £12 [£11] denotes the 3x3 intersite [on-site] self- 
energy matrix in the space of the ti g orbitals. This is then 
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combined with i?i^;(k) in order to obtain the Green's 
function at a given k-point. After summation over k, 
the intra-dimer block of the Green's function G^j, is 
extracted and used in the C-DMFT self-consistency con- 
dition. The 6-orbital impurity problem is solved by a 
numerically exact quantum Monte Carlo (QMC) scheme 
using up to 100 slices in imaginary time at tempera- 
tures down to 770 K. From the Green's function on the 
imaginary-time axis we calculate the spectral function by 
using the maximum entropy technique. 

We first discuss our results for the rutile phase (Fig.QJ. 
In this phase, we found that the results of single-site and 
cluster-DMFT calculations are indistinguishable for all 
practical purposes. Correlations reduce the total (occu- 
pied and unoccupied) bandwidth corresponding to the 
coherent part of the spectral density from 2.59 eV in 
LDA down to about 1.8 eV (see Fig. This is consis- 
tent with the quasiparticle weight that we extract from 
the self-energy: Z ~ 0.66. Hence, the metallic phase of 
VO2 can be characterized as a metal with an intermediate 
level of correlations. We observe however, that the occu- 
pied bandwidth is barely modified (~ 0.5eV in both the 
LDA and our results). This finding is in agreement with 
photoemission experiments [IRIl7| |. and our work demon- 
strates that it is compatible with a local self-energy. A 
prominent quasi-particle coherence peak is found close 
to the Fermi level. This has recently been demonstrated 
experimentally in high photon energy bulk-sensitive pho- 
toemission experiments 17] (and is also suggested by 
low-photon ener gy PES provided surface contributions 
are substracted |lj|). Hubbard bands are apparent at 
higher energies, with a rather weak lower Hubbard band 
about —1.5 eV below the Fermi level and a more pro- 
nounced upper Hubbard band at about +2.5 — 3 eV. 
The former has been observed in photoemission exper- 
iments 0, 0, 0] , while the latter constitutes a predic- 
tion for BIS experiments. Finally, correlations result in 
a slight increase of the occupancy of the a\ g band (0.42) 
relative to the ones (0.29, 0.29), in comparison to LDA: 
(0.36; 0.32, 0.32). 

Figure [2] displays the spectral functions for the mono- 
clinic phase calculated within cluster-DMFT in compar- 
ison with the LDA DOS. The key point is that inclusion 
of non-local self-energy effects succeeds in opening up a 
gap of about 0.6 eV, in reasonable agreement with exper- 
iments. Our calculations also yield a large redistribution 
of the electronic occupancies in favor of the a\ g orbital 
which now carries 0.8 electrons per vanadium (with only 
~ (0.1,0.1) remaining in the e g orbitals). This charge 
redistribution is a common feature of many models of 
VO2, whether correlations- or band-driven, and is also 
observed in experiments. In the present context, the fact 
that the single electron occupies almost entirely the a\ g 
orbital, together with the fact that U is larger than the 
intra-dimer hopping (itself much bigger than other hop- 
pings), means that the ground-state of each dimer is close 
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FIG. 1: Spectral function for the rutile phase as calculated 
within DMFT with U = 4eV, J = 0.68eV (solid lines) in 
comparison to the LDA DOS (dashed lines). The red (blue) 
lines show the partial contributions of the ai s (ej) bands. 
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FIG. 2: Spectral function for the monoclinic phase as calcu- 
lated within cluster-DMFT with U = 4eV, J = 0.68eV (solid 
lines) in comparison to the LDA DOS (dashed lines). The 
red [blue] lines show the partial contributions of the ai s [e^ 
(averaged)] bands. 



to the Heitler-London limit. In this limit, one has two 
electrons with opposite spins on each site forming a sin- 
glet state, rather than an uncorrelated wave-function in 
which two electrons (per dimer) are placed in the bond- 
ing combination of atomic orbitals (hence a large double 
occupancy). The transition into the insulating state is fa- 
cilitated by the Heitler-London stabilisation energy. This 
confirms the early proposal of Sommers and Doniach |j| . 
As a result, the dominant excitation when adding an elec- 
tron in the a\ g orbital costs an energy which is set by U. 
This is apparent on our spectra: the antibonding combi- 
nation of a\ g orbitals corresponding to the narrow peak 
at ~ +1.5 eV present in the LDA calculation is replaced 
by a broad upper Hubbard band centered at ~ +2.2 eV 
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(the precise value of course depends on the choice of U). 
Correspondingly, there is a (weak) lower Hubbard band 
at ~ —1.8 eV. The low-energy nature of this singlet insu- 
lator, however, is quite different from that of a standard 
Mott insulator in which local moments are formed in the 
insulating state. This is manifested by the fact that the 
low-frequency behaviour of the on-site component of the 
self-energy associated with the a\ g orbital (Fig. EJ) is lin- 
ear in frequency En = £n(0) + (1 — l/Zi)u> + • • • , in 
contrast to the local-moment Mott insulator in which 
En l/^- 1 diverges. As a result, the spectrum dis- 
plays a narrow a\ g peak, at an energy ~ —0.8 eV below 
the gap, which carries most of the spectral weight for 
u> < 0. This peak should not be interpreted as an inco- 
herent lower Hubbard band, but rather as a quasiparti- 
cle which has been gapped out (Zi can be interpreted as 
the spectral weight of the gapped low-energy quasiparti- 
cle in the insulator). Hence, at low-energy, the physics 
is that of a renormalised Peierls insulator (analogously, 
a correlated Kondo insulator can be viewed as a renor- 
malised hybridisation-gap insulator at low-energy). The 
bonding-antibonding splitting is renormalized down by 
correlations. Indeed, a weaker a\ g peak is visible in our 
spectra at the upper gap edge, at an energy considerably 
smaller than the antibonding peak in the LDA DOS. The 
e g band, in contrast to the a\ g , is weakly correlated, as 
evident from the self-energy in Fig. and expected from 
the low electron occupancy. Its bottom lies in the same 
energy range than the "renormalised" antibonding peak, 
so that the gap can as well be considered to open between 
the aig and the band. 
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FIG. 3: Self-energies £(iw„) from the C-DMFT calculation 
(for U=4 eV, J=0.68 eV) for the Mi phase, at low (imag- 
inary) frequency. Circles [triangles]: imaginary part of the 
on-site diagonal a\ g [e^] component. Diamonds [stars]: real 
[imaginary] part of the inter-site ai fl component. 

Our findings for the spectral function explain the re- 
cent photoemission data of T. Koethe et al. [ijj. These 
authors made the puzzling observation that a prominent 
peak below the gap appears in the insulator but at an 



energy (~ —0.8 eV) which is shifted towards positive en- 
ergy in comparison to the (weak) lower Hubbard band 
measured in the metallic phase (~ —1.2 eV). We propose 
that this peak is actually the coherent quasiparticle peak 
at the bottom of the gap, rather than the lower Hubbard 
band. The latter is present at higher (negative) energy 
but remains rather weak even in the insulator. Fig. 0] 
compares the total t2g spectral functions calculated for 
the metallic and the insulating phases. 
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FIG. 4: Spectral function of the t2 S orbitals for rutile and 
monoclinic phases as calculated within DMFT/cluster-DMFT 
with U = 4eV, J = 0.68eV. The Hubbard band in the rutile 
phase is located at about —1.2 eV, whereas the pronounced 
peak at — 0.8eV in the monoclinic phase corresponds to the 
gapped quasi-particle (see text). 

Finally, let us mention that we have discussed the re- 
sults obtained for a specific choice of the interaction pa- 
rameters U = 4 eV and J = 0.68 eV. We have actually 
performed calculations for other choices. In particular, 
we found that it is possible to stabilize an insulating state 
within C-DMFT for smaller values of U (e.g U = 2 eV), 
but only if the Hund's coupling J is taken to be small 
(in which case it is easier to redistribute the charge to- 
wards the a\g orbital). Also, we emphasize that a single- 
site DMFT calculation can lead to an insulating state for 
rather large values of U (U > 5 eV). However, this in- 
sulating solution is a conventional Mott insulator with a 
local moment, and therefore would display a large mag- 
netic susceptibility which does not correspond to the ac- 
tual physics of insulating V02- 

In conclusion, we have presented LDA+CDMFT cal- 
culations for vanadium dioxide. Both the metallic rutile 
phase and the insulating monoclinic phase are correctly 
captured by this approach. 
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